


Institutional Archive of the Naval Postgraduate School 


Calhoun: The NPS Institutional Archive 
DSpace Repository 


Theses and Dissertations l. Thesis and Dissertation Collection, all items 


1985 


The effectiveness of heat exchangers with 
one shell pass and three tube passes. 


O'Hare, Mark S. 


http://hdl.handle.net/10945/21289 


Downloaded from NPS Archive: Calhoun 


| Calhoun is the Naval Postgraduate School's public access digital repository for 
D U DLEY research materials and institutional publications created by the NPS community. 
22 Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS's first 
KNOX appointed — and published — scholarly author. 





LIBRARY Dudley Knox Library / Naval Postgraduate School 
411 Dyer Road / 1 University Circle 


http://www.nps.edu/library Monterey, California USA 93943 










































































Ал . 72.4, 774 ЛЬ, Р ЦА 
гс"... є А E VN LY TEQUE 
ата" 3 АД ЧОГО ГУ | у eh tt АДА 
И oa АЛ Ка. Е Rum Ў КАЕ АК ТО ek vy o pou) 
Ер" 4, л КИЛАТ КАЛА Р а ИНА ДАВА САЗА КАСА а бү ЭЛ. Y Заах 
сл 220. Ор ре пин лан а Ско ха сти | ОА ЭМ 
217 M" Ы в. SAM rM RS E кұ ах, пам a Ау: Ton М. 
V^ 214: м " и а аа [NW] л. № а + ~ ч А 524, АА 91 f те » > 
+ ы La + А АЈС У уочи 44 а Ши M 
Ћ, ээ А: 4% 3 2” ҮЭ " т Ч We КВА *9- 9 t х ^ 
р 1 H 1 В “3 и - 
А7 ONSE NE eh ET, ХУ UG) DOES Wn, 
ГА - м y ^ П А 14, A ы whea А. 7 я көз 
I А TAI АЧГ [К ^ (ОЧ A PA. А. 
ы “ 1 у LI LI 1 PN & ч TV ШУУ Ae Д. УБ | 
5 2 Y З МИ ДА a t è N NA. ' 545 е) * A "ҮҮ, 
б E Ч . B мг» 47 А» 
ЖЭЛ Й RN TES oU RUP S Mr eh n m e Ааа 
' y Ци 5 2 Tu гү үлгэр h ғ 4 
" a 8 М к 5 b pr 4% DM АА, Y А vy Ree % М; ъ 
? А А П БУ Ми о", ~ ГУЗ "WO A М 4 
зон D 1 с^ КТ УГ ч эт PUOI EC S КАЖА КК 
TT . 3) Y VS A HT Ж PT LEY DUM. n "P 
мі 4 | . Ц к; ^ “ м “ Ал ru p TAA я 
T де » 4 *) "84, x М Ч Эгч Де АД ТОК 
( + қ 1 FS жас ДАДАТ АА ча CUM PUR Де ВАМ 
^ Uu И“ d Ж т ТҰСЫ NUN prac hy QUU Nee A AUC EM ESAE 
, t е D» b v», aca - ТҰ (Үг cm А | Оо S. oen л, А ҮТҮЕ; 
Ч ус о АД zo rua S mA. Ys A CCS Ar AEn eV aT ч 
, А А. И т. ЖА” АЧ < OU wd M eat: Миа TA У Нара МА 
4 МР у AAA Soe el ey el "Twy сү AD КОСУ О 
D , 1 ТЭЭГ: ЭГ ТҮ г і „а 2 ; „Л, (Ж IY 
Ж; " К 3 t 44 А Ж Ка | П қ Ұқ B ЭХ, да ааа EN 2. ы) POS Да уан ЁС 
Y В лы? “уы” г дайы" 5577 р 5 ұң Жа 
» | ы * UN. Me С , a > YT V: Ка 4 чак КМ М КЕ Be Ses NAT а ТА ќи jd 
П Е L EXPL жэл үр ud ҮЛЭГ ЈА ОТ л Vic Pong RS A AUTEM ын ево А E 
< Ж г ? я Ха $ 44 $ + 7 [ еу CUN Pa Ч н tF R. dip ДО: а = ДА АЧ 
E қ У | 4 ч а Ад, а: м Мы? ЗАК i ДС НА; ALL рут ау тора ае и рам Но, 
- Yt. сүү ДЕКА сно а бо o) DN Mr Ae Ri КАСА A IQQ К О Vela ^ , Колева . 
rw П | DES 8 р Ша. лаз | | | & pu j Ч 
uve Loi а „ыйы АЕ А СУ ТОА уа СИА ИУ LE d D em А иг. А ren 
Ж П " ` АЖА Үт. %%4/% 4» мр NN VP г. г м Эл, Mar au Pv Хана п » t 1 | 
- & v , Г t » ХҮХЭЭЛ Ч р E. " МОДА КОЖИ а, Яне, Worry e Y aa LATA RÀ. ы, ЗА Ал 
з % d ы Б PW AX Имам шэг Кра ми Ай Ya es Й b 3 үл НАИ $ п» tent Ae Sy] 1 АА, Цан р 79 : 
3 ЖОЛОО Е қалала UE Seo Соба шэн ЛА ЫР ДӘ 
" А А и ELM МА АТР АЛЛАА Ср, АМАЛ ҮЛ ГЭ БУЙЛ i 3, Када Ах УЗ а че ДТ уча UN MA Л сан 
: ЛЭ < ч <. n P ay А МАЗ хаха Lu SC А Мк» х УХАА", Mati! v №; Чэ Я А Ч: AUS mA vd Tw Tp 
М ДА wey мал ОК Л А RC » UN res LOC SOUS ae ty 
А Н пес А Oe a er я TAM i Гос а ыы ! AIC Oy" SA ёл ва 9 эг hy, Hel te, Be, р 
г ' Эг A пача SI конот ад Ұқ VAT A S es А: 1 + Негри ^ 
x D 4 ЛАЛА ГТ T ur КТ ah € yp ats "e ~ Pe DUO оту ма Јо 4 ЗЭ SRSA е А) 
y [1 Р « E LS ip (і и 1 5% га, и V m Де м 9”, EXIST ОРЧ M нэгд x "у « 25% Түлек 4 
2 d + 29 5.94 №, > Varian} SPR D nd nr Г 
Г 27 ЗА И uu НГ ЈЕ ем SVS arte НАТА ДКА, 4 МА Xe] НИС м 
н Ы f DE Y “е АРЦ EC M АЖИН ДОЛ BOO COLS AREA SANA AR J 
р ХАЛ" ак 4 Мм р | 
. 3 t «X l^ zn ч x , e E ҰМТ; be t АК Л “ 2); X М ы; Ка AA e к 
, А М АЖА МА ААА ае TUA v 
„2 E SA тол Ул ми ЗУ 54 С „ку БН УСКЕ | М АНА Из ын РУТ 4 7 
и ч қ A ШАН Ај Vets E Ээ МА Ж. DAT! UO rh PIE LAUS , МИНУ 
мат " 127% САЛЖ, C PU AX T ро ке КА ҚАҒА сауу НИЈЕ А 
П " ... ) 35 11" ГҮЙ” ГА хх қал” 5 | JA 9 та ЗАС An ен Li i iat A Ur Ги 3 уян З, г " 
А А 27 D Br , Ч! R ОЛ, 4: аи 4 г. d Ч 4 ix AM "3 ИЗА Dy АЖ (39 M | 7 
| % р. - AL 8 де“ Ми ^? и 
A Eun UL WEE Y Eun Wh nr ” HS | S Ud C ој Tue зму 
У м 1 У з 4 ь ИХ, № d чо EN 25 " 74: 19 ныг 76: 1 АР АН g ы 4 Yad УАН te АД 
ls у | 1 LUAM DCN ЭЛ т A Wh oh? пад, а Мөст М 4 РАЈ ^ y 
ГІ; MS 3 M. Ly 1 К 1 ALPEN" ғал LA A LP 4% “фр Ха n 2% а 479: С + й Ке de До 
| Водич ун ^ қ рани Еч 
ман, 2 од Й ТО УУ ДААХ ЛЭН СКЕЧ от ОҚЫР Аға 
HA ка Ж . : Y MEA EV LULA o DO Киа НЕЛИНА 4 
UJ ^ ONE ? - З ла EDU TEE qa У А. қ Ч, qd 3:9 dv. ) 29996 la. КҮ PINE VIN DP жалы 
! Ц ^ E О 47% О 4 қ жы 2, 4 ET. 1 17 E : AT Ре > Е f> ^ “ч ХА 
А А Ut ~ " EC РО "T 7 1 У "^ Р ! 
3 За. № è Y Е П 
$ цагг + | и ча - 
x , 
1 LM - 


(3 ч Hi А e 
V PAAR aM MU, ААХ Бад хаҡ лах Ў 
лу", МА АА қ с | г ұз КАБА аа MIX. > ; x КЧ ОНОХ т БЕЧА 
бум? n АКТА: Ані $ P 2%4 ДАШКА А У. д Мақ, гү б Хы ои г к P Ж м КУ У 
Е 1 учено лу ни ца пок АА UE LR NER AR КАА: МИ РА А У ЕКА, а дра 
[] “| g ч Ы M i e 117121 Др 119 " 8 ы / г ^ | ^ 2 Ы a им ы 4 7 с 
Ж Й va. ма! а 4 PAPE A му. ЗЭ. D LA ^ a 6 ка ; МА dM Shi iN ob A CAN цал 2! 
^ 1 5 МП Р 1797 
2 CEDE. LEA А UA ДЕКА таа ^ 


ЕУ 
к: 
M 5 Je ГАЯ А, 









































ү” 
я й Xt MA > 
КЕ ЗАМ НУ ЩЕ Ска Sage А Ера Tct, A 
4. ад. ЫЗ, Ји ТЗ У А k АДА Хо b “~ ра TU 22 
J З ay ча ч i 4% мэ we d Ы: A b. р M 
UP , Ч Ха АА Ads "X M dat Ar lA e Dra D ?, RASA A. 125170 7 
3 % ^ D 2 ЈУ ^ "t ~" г. «гай yu қы” ~, - 
А П 3 Aa ys мэтэ А ХОГ” AA a A oa 4,5; te e je a, kee 
121 үл СО Ou PN. 77% APR ~ л C tur Vs 
5 4, ЖЭ, ГОЛДОГ ЭЛ А" ДО ЕН Эг PE Y Pr ACA ak a MCN А S 
"LU ч. ЦЭГ" Aras Кр рее, M | Көлі ұз 
И | bl: PE * À < пі аб а а ва А И А 
3 LJ % ~ Ка Mx M 12-9 ПИЦА я af 7 мэд 
è п A 4 иа еи ОДА ш, ЧАГА ДД, ХА УУ: ЭГЧ: 
Ёс ТЕЛ " , 9 ГА АА EAT Yo P A Рим 
p 7 Е, BE T. Ж Н 4 A YX ми AP 45 УК iiA її, ^ по ule 
| d ' - АЕ 1 АИЫ CP HEX Os 4% 5 + РАМА C ad 
4 1 ч ОЙ А kYYw А? ЖАГ tQ Маг су ка 
АХУЛ А СЕ АИИ ELO и Са арии РАНАМА: ААКУ УЕ 
П Ж” Со ЧИ АА, гор A та У 53! салад 
1 | тэр ко CIAO EUN SOUL reba hi: 
Й ак М Lr EROS UA “уз чл Рк ЗИ ЖКС НЕТ 44 
| “ы Элл СЪД м Ч Ға"! У 
Й Г з ' Wa | 
L] D 47 Е 
Г : Ы % ч $c к 
Ж ` A ҮЕ, a 
ЭР P nf Эг и y^ 
“% > м” - 9) ^ М А 
, LI y h iX му * 
О % г f Ц E m ЖАТА е Р ГУР; қ 
РУЧЕК ДЕ es АСА бр ДОК 
Ка қым” EC CS Т 
Y SUP ч 
" А, 
г 
А 27 А X | 
M , 7. 
Га я E 
Н О 
[ 
b 
2 
5 Г. А Ds 
- ла 13: 
E П 21) 
E 
- ото 
i 
A 
НО Ги МА 202 
"TES n У ^ Е и "Ж; ія іч 
УУ ; У РИА АЕ р 
Й $ As ROS) ТАЛАСТА АРУ at EREL И TE E QU e M P C 
> ТАС: у 2-3 bord ! ИДЕ Хак. CEU ИХ PA АЧ 
і - EP J б af PNE, a v? б,” ДУ, EXE INS “үс ЗА aD C DORT B7. 2% 7 
М У - 25 PA РИ 1 CN ib LY 4 SAG НЕ V. x ys es 
СЫ Еее НИТ КЛ АД С АЯ ка М РА УЧ RAS PTS. 
| n 2 А Vul алы uet a АЕ. 1 Mae Sá хе 
1.:/% ГАЙ 5 РА КА dA EM à 1 #3 РМ та їх 2 
А Џ * + а 4 ү ЕУ be | Ул (2 
LEUTEN ӨР ПРЕ и ҚЫҚ Дай 
f ovd 
Li И з 
POR А 
fa 
[ г 
» “ч Га 
Н г 
T 
4 
„Бад 
y 
5.4 | 
i , 
2 Ф 
ка. 
4 
AUAM 
2 
ЖАТ П 
Ls AR] 
TA, 
E d F 
Ц 
"n 
г 
LI 
' ТЫШ; 
é 
% 
? 
e f 
LI 
M 4 
г 
4 1 Papà 
Фи ХАД биа Y 
A ч ү SU 4 ^ 
Ӯ 4 Й 8017 РИ 
» КАТА ра J 
ж 3 Ur LIU d 
НАЈ | 
й XA o A OE 
| И р 
^ Р 3 : кр ја 4 р X УРА А 
у i б ААА Ө, x У. кр й а: А Т КАН е М 
Р y Е s қ NUN AME X fautes Sosa aii я A AUI пи БТА е И 3 Ма 12 
ft yl LUE y a У td Ando Fu EC За ЗАС наи Ан ЖЫ УА a и +: prag г РЕ нх 
зы 4 | 3%» 74 , ќе ААА, P Ж LRL J fy А ; БҮЛ, и 1 МЕ т ғ мэ Фет Ын Te Le 
и Y MT p of FF ы і d JN D у (16) 1 4” M 7:6. АЈ үл чи" m 
M НО СО А А | S ae tS DOE CO TUR eee 
1 д É "ne t s ў Ў 1 ит У | БЕКЕТІ COD. J 
МА ДЕ Е TAN х | ; | И > t PAP Ми зам АДА б 
| т 5 Пе Аа бо AMETE e У 4 ОНЫ 
К JL Lb А .. ко АЈ Гел ^ AM. Tr T В И Bos A ч Ч P UA +" ч 225 АА 
: ДАНА У pun eene Дада | TR Upper SE 
| me ка ети Рачу." ди ШУА Ц емч ружи лон Ын дэ ЦА таз: M п IOS AR 
й ^s. M I E РЕД Ааа 5 Я ш, АА а ақа И РСС а ВАТ РАД erret babe e 3. А je d ч, po 
: EUR АО ЖА, МЕ d ^ МА E + >» | Ч оа EEN, д мэ + 3 а, я 
d р i 1. р : > 1 я RA з зи НА. д D 1: ЧИН (Р 14 У 4 Ф м б у "р 27:40 са ф MN 
5 А етеу АК, ~ (1 Сү D Aci гү? E X n А У 
1 E] ji g^ ”% 
. П б ~ (3 DE r ы Y, dr Кај s dV 5 
Зе i (0 n pL О “ бу ГЖ ЛА, D $ 
7 қ ЯС 1 Ч 
» 4 р б, 
| S t 
Н 1) 
Л 
П "У 
' » 
үй 0 
4 
/ a . 
LI | 
г р а 
Г. 
А " 






814 р H E 2 ЗУ 

7 К Ч в UM про - M 1 Ч Б [A Ж L s 

І ; 4 AAT REL Tn оче) Ма А А МА AUTE ТАЛ ЭХ А о ^ БА 5, 
ња. үү КОЛАТА ДА Уард қ А му > "€; 
BONN ЛИПИ ЖЭ) 
я " Rm T 3X л А 533 EC Ae ОЙ. " 
АА omy | 17 ТЕ bat a te и МЗТ, ава 
5 : (32 4 . 





DUDLEY KO T LJU # 
МАМА, РОГ | 
МОНТЁХСҮ, Cide. 








NAVAL POSTGRADUATE SCHOOL 


Monterey, Galifornia 





THE SJS 


CHE TEFFECTIVENESS OF HEAT EXCHANGERS 
WITH ONE SHELL PASS AND 
THREE TUBETPASSES 


by 


Mark S. O'Hare 


June 1985 


Шнесіс 40150016 Allan D. Kraus 





Approved for public release; distribution is unlimited 





SECURITY CLASSIFICATION OF THIS PAGE (When Data Entered) 


REPORT DOCUMENTATION PAGE BER EE UE SE ам 


1. REPORT NUMBER 2. GOVT ACCESSION МО. 3. RECIPIENT'S CATALOG NUMBER 


4. TITLE (and Subtitle) 5. TYPE OF REPORT & PERIOD COVERED 
The Effectiveness of Heat Exchangers Master's Thesis; 
With One Shell Pass and Three Tube June 1985 


Passes 6. PERFORMING ORG. REPORT NUMBER 


AUTHOR(s) 8. CONTRACT OR GRANT NUMBER(s) 


















Mark S. O'Hare 


. PERFORMING ORGANIZATION NAME AND ADDRESS 10. PROGRAM ELEMENT, PROJECT, TASK 


Naval Postgraduate School 
Monterey, California 93943-5100 


AREA & WORK UNIT NUMBERS 


11. CONTROLLING OFFICE NAME AND ADDRESS 12. REPORT ОАТЕ 
June 1985 
Naval Postgraduate School ИЕ ВЕНЕ ТТ” 
Monterey, California 93943-5100 | 204 


14. MONITORING AGENCY NAME & ADDRESS (If different from Controlling Office) 15. SECURITY CLASS. (of this report) | 
Naval Postgraduate School | Unclassified | 


Monterey, California 93943-5100 
15а. Е AT ON DOWNGRADING 


ағасы 9 


16. DISTRIBUTION STATEMENT (of this Report) 


Approved for public release; distribution is unlimited 


117. DISTRIBUTION STATEMENT (of the abstract entered in Block 20, If different from Report) 


18. SUPPLEMENTARY NOTES 


19. KEY WORDS (Continue on reverse side if necessary and identify by block number) 


Heat Exchanger 
Effectiveness N М Method 
Three Tube Passes 


20. ABSTRACT (Continue on reverse side If necessary and identify by bfock number) 


Heat exchangers with one shell pass and n tube passes are 
Often referred to as l-n exchangers. The heat transfer 
literature contains many references to studies of l-n exchangers 
when n is even but apparently little work has been done with 
respect to the l-n exchanger when n is odd. This thesis 
greatly expands the theoretical study of l-n exchangers with 
n being odd. While a completely closed form solution was 


DD , FORM 1473 DITION OF | NOV 65 1S OBSOLETE 


JAN 73 
ео аш SECURITY CLASSIFICATION OF THIS PAGE (When Data Entered) 


SECURITY CLASSIFICATION OF THIS PAGE (When Data Entered) 


207 (Continued) 


found to be unfeasible, a polynomial approximation has 

been developed that yields the effectiveness (=) of the 

two possible arrangements of the 1-3 exchanger as a 
function of the capacity rate ratio (R) and the number of 
transfer units (N E It is also shown that the effective- 


ness of the arrangement with two counterflow and one 
parallel flow tube side passes exceeds that of some of the 


l-n exchangers with n even. 


S/N 0102- LF- 014- 6601 


—  Ó—— — à Дана ња ани 
2 SECURITY CLASSIFICATION OF THIS PAGE(When Data Entered) 


Approved for public release; distribution is unlimited 


The Effectiveness of Heat Exchangers 
With One Shell Pass and 
Three Tube Passes 


hy 


Mark S. O'Hare 
Lieutenant Commander, United States Navy 
B.S., United States Naval Academy, 1976 


Submitted in partial fulfillment of the 
requirements for the degree of 


MASTER OF SCIENCE IN MECHANICAL ENGINEERING 
from the 


NAVAL POSTGRADUATE SCHOOL 
June 1985 


ABSTRACT 


Heat exchangers with one shell pass and n tube passes 
are often referred to as 1-п exchangers. The heat transfer 
literature contains many references to studies of 1-п 
exchangers when n is even but apparently little work has 
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This thesis greatly expands the theoretical study of 1-п 
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mation has been developed that yields the effectiveness (ec) 
of the two possible arrangements of the 1-3 exchanger as a 
function of the capacity rate ratio (R) and the number of 
transfer units (Neg). It is also shown that the effective- 
ness of the arrangement with two counterflow and one 
parallel flow tube side passes exceeds that of some of the 
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NOMENCLATURE 


English Letter Symbols 


A 
Am 
Ao 


А1 


Ад 


Аз 


А4 


А5 


Exchanger heat-transfer surface, sq m 
Coefficient of the mth value dimensionless 
Coefficient, dimensionless 


15$ order coefficient to be multiplied Бу Му, 
dimensionless 


2nd order coefficient to be multiplied by МЕ, 
dimensionless 


3rd order coefficient to be multiplied by М№+13, 
dimensionless 


4th order coefficient to be multiplied by МЕ 
dimensionless 


oth order coefficient to be multiplied by Ney”, 
dimensionless 


Exchanger heat~-transfer surface, sq m/m 


Capacity rate, W°/K. Also designates dimensionless 
arbitrary constant 


Specific heat at constant pressure of cold fluid, 
J/kg?K 


Specific heat at constant pressure of hot fluid, 
J/kg?K 


Empirical value of effectiveness (computer generated), 
dimensionless 


Error. Also used as the exponential function 


Logarithmic mean temperature difference correction 
factor, dimensionless 


Exchanger length, m 
Number of effectiveness empirical data points used to 


determine a curve for R, dimensionless 


120 


Number of tube passes, dimensionless. Also number of 
equations, dimensionless 


A related Ney per unit length hot side, m71 

A related Nty per unit length cold side, т-1 
Number of transfer units, dimensionless 
Temperature group, dimensionless 

Total rate of heat transfer, W 

Maximum total rate of heat transfer, W 

Capacity rate ratio, dimensionless 

Temperature group, dimensionless 

Sum of the squares of the residuals, dimensionless 
Hot fluid temperature, ?C 

Particular integral, dimensionless 

Hot fluid temperature in, ?С 

Hot fluid temperature out, “С 

Cold fluid temperature in, ?С 

Cold fluid temperature out, ?C 

Cold fluid temperature list pass, °C 

Cold fluid temperature between 1st and 2nd passes 
Cold fluid temperature 2nd pass, °C 

Cold fluid temperature between 2nd and 3rd passes 
Cold fluid temperature 3rd pass, “С 

Overall heat transfer coefficient, W/m? - ?С 


Mass flow, kg/sec. Also the product of wand L, 
dimensionless 


length coordinate, m. Also used to represent 
a constant value in a Sequence 


Ша 


y = Sum of mth degree polynomial, defined by eq. (53), 
dimensionless 

7 = A product of z and L, dimensionless 

2, = A related Ney per unit length, hot side, 1/m 

Greek Letter Symbols 

a = Root of auxiliary differential equation, 1/m 

Є = Exchanger effectiveness, dimensionless 

À = Combination of variables defined by equation (11), 
dimensionless 

ш = A related Nty per unit length, hot side, 1/m 

ф = А combination of terms defined by eq. (38), 
dimensionless 

> = Summation, dimensionless 

c? = Variance, dimensionless 

да = Mean temperature difference for exchanger, °C 

д = Indicates partial derivative, dimensionless 

Subscripts 

С = Со1а fluid 

ћ = Hot fluid 

i,j,k = Values in a sequence 

m = Degree or order, an exponent 

| = inlet 

2 = outlet 


Special Symbols 


[A] = 


An m x n matrix, symmetric 


Ез 


[K] 


[L] 


An mx n matrix, symmetric 


A lower triangular matrix 
[Т] = Ап m x 1 vector 


[0] 


A null vector 
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Io INTRODUCETON 


A. BACKGROUND 

When analyzing the standard counterflow heat exchanger, 
it becomes apparent that, from a practical standpoint, it is 
often difficult to obtain a high velocity for one of the 
fluids when this fluid is constrained to flow through all of 
the tubes in a single pass. This leads to a possibility of 
a low overall heat transfer coefficient which cancels the 
advantage of the high logarithmic mean temperature 
difference which is obtainable in true counterflow. 

The quest for flow arrangements for increased heat 
recovery has led to arrangements that yield increased tube- 
side velocities and higher overall heat transfer 
coefficients even at the expense of a departure from the 
ideal true counterflow arrangement. Thus, the design may be 
modified so that the tube side fluid is carried through 
fractions of the tubes consecutively. 

Heat exchangers of this type with one shell pass and n 
tube passes are often referred to as i-n exchangers. These 
exchangers, such as the one shell pass, two-tube pass (1-2) 
parallel-counterflow exchanger (see Figure 1.1), are 
configured such that all of the tube side fluid flows 
through the two halves of the the tubes successively. A 


Single channel is employed with a partition to permit the 
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Тнотин (T4) 
PARTITION 






Ссоголм к) 


{согр,ошт (tj) 
FLUID RETURN 


Тнот,очт (T3) 


Figure 1.1 1-2 Parallel-Counterflow Exchanger 


entry and exit of the tube side fluid from the same channel. 
Note the baffles used to induce turbulence causing the 
liquid to flow through the shell at right angles to the axes 
of the tubes thus helping to create a higher shell side 
velocity with higher shell side heat-transfer coefficients. 

To date, much work has been done on finding the true 
logarithmic mean temperature difference for heat exchangers 
with an even number of tube passes. Little work, however, 
has been done with regard to heat exchangers having an odd 
number of passes. This is primarily due to the method 
employed in deriving an analytical solution to measure the 
overall effectiveness of a heat exchanger. 

Exhangers with an even number of tube passes often 
present a configuration problem especially in a marine appli- 


cation where the inlet and outlet of the cooling fluid must 


ШЕ) 


be one the same side of the exchanger header (see Figure 
1.1). This particular problem could be alleviated by going 
to heat exchangers with an odd number of tube passes. 
Currently, precise mathematical expressions for the 
effectiveness of the 1-3 exchangers do not exist. Hence, a 
theoretical examination is reported on here which considers 
the effectiveness of the 1-3 parallel-counter flow exchanger 
which is shown in Figure 1.2. 

Before doing this it is important to mention the work 
that lead to the effectiveness method and the development of 


work on heat exchangers with an odd number of tube passes. 


нотам (т; 
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Figure 1.2 1-3 One Shell Pass Three Tube Pass Exchanger 
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This will be considered in Section II. Кегп [Ref. 1: 

pp. 224-226], makes the interesting point that the optimum 
exchanger requires an exchanger capable of providing the 
optimum fluid-flow velocities on the shell as well as the 
tube sides. This might frequently entail the use of an odd 


number of tube passes or an odd tube length. 


B. WHY EFFECTIVENESS AS A FUNCTION OF Ney 
It is also noted that Kays and London [Ref. 2: pp. 
24-29] indicated that the effectiveness as a function of 
Мец Се - Ntu) method is the favored approach for evaluating 
a heat exchanger's performance because: 
1. The effectiveness value stands alone as a dependent 
variable and should not appear directly in the 


abscissa and indirectly in the ordinate of a graphical 
display. 


№ 


The log-mean difference equation misleadingly 
simplifies the notion of what is involved in heat 
exchanger design theory, since the implication is that 
only a rate equation is required. 


3. Тһе e - Neq, approach simplifies the algebra involved 
in predicting the performance of complex flow 
arrangements. 


4. The more meaningful arguments are related to ease of 
use in design work. Two prime examples of these are: 


a) Given the overal heat transfer coefficient, U, the 
two fluid capacity rates, Co and Сү, апа the 
terminal temperatures, determine the required 
surface area, A. 


b) С1уеп А, П, Се, Су and the inlet temperatures of 
both streams, determine the outlet temperatures. 
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II. THE DEVELOPMENT OF THE EFFECTIVENESS METHOD 


A. LITERATURE SURVEY 

It was Nagle [Ref. 3: pp. 604-609], in 1931, who 
credited Davis [Ref. 4] with a simplified method for 
computing actual temperature differences between two heat- 
interchanging streams which depart from true counter or 
concurrent (parallel) flow. This is now the familiar "F 
factor" method which expresses the actual mean temperature 
difference ðm in q = UAðm as a fraction F of the counterflow 
logarithmic mean temperature difference, LTMD, фе via 6m = 
Е Өле, 

The example of initial interest was the 1-2 exchanger 
with a single shell pass and two continuous tube passes in 
counter and concurrent flow with it. The method involved 
derivation of the actual temperature difference for the flow 
pattern and formed the ratio F = 6)/ һе: This familiar LMTD 
correction factor was plotted conveniently as functions of 
the effectiveness, є, and the capacity rate ratio R with R 
as a parameter. These mean temperature difference 
correction charts are available for many flow arrangements 
[Ref. 1: pp. 829-833 and Ref. 5]. The effectiveness, e 
(often called P or S), is always the cold fluid 
effectiveness and R is always the capacity rate ratio of 
cold fluid to hot ЕТО 
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Nagle detailed assumptions and derivations for the 1-2, 
1-4 and 1-6 exchangers. The F factors were obtained by 
Nagle through graphical integration and were accompanied by 
the comment that F factors for the 1-2 exchanger could be 
applied with negligible error to 1-4 and 1-6 exchangers. 
Underwood [Ref. 6: pp. 145-148] rederived the equations of 
Nagle for 1-2 and 1-4 exchangers to eliminate the need for 
obtaining F factors by graphical integration. 

Bowman [Ref. 7: pp. 541-544 ] pointed out that for a 
very large or infinite number of tube passes, the F factor 
approached, as a limit, its value in crossflow with both 
fluids completely mixed. ТЕ was further stated that even at 
the limit, the F factors were only 1 to 2 percent lower than 
those for the 1-2 exchanger. A previous paper by Kraus and 
Kern [Ref. 8] did not confirm the generalization that 1-п 
exchangers differed only negligibly from the 1-2 exchanger 
although this lack of confirmation was obtained on an 
= = Е(В, №.) basis. Moreover, the Kraus-Kern work does not 
confirm the generalizations on an F = f(R, Ney, £) basis. 

From the standpoint of usefulness and good accuracy, it 
is essential that F factors, if they are to be used in 
preference to в + f(R, Ntu, flow arrangement), бе obtained 
with precision. Plots of F = f(R, e, = Por S) [Ref. 1: 
pp. 829-833 and Ref. 5] show that the curves for particular 
values of R approach infinite slope as F decreases. While 


this can be partially alleviated by restricting R < 1.0 (a 
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constraint used in the e« = f(R, Ney, flow arrangement) 
approach), it is seen that small errors in the interpolation 
for R or e = P or S can result in large fluctuations in the 
value of F. 

In a comprehensive paper, Bowman, Mueller and Nagle 
[Ref. 9: pp. 283-294] presented graphs of F factors for 
shells with one through six shell passes and numbers of 
continuous tube passes respectively double the number of 
shell passes. In view of the earlier references to Nagle 
and Bowman, it should be noted that F factors were computed 
for the 1-2 exchanger in [Ref. 9: pp. 283-294] using the 
equations of Underwood [Ref. 6: pp. 145-158]. 

Ten Broeck [Ref. 10: pp. 1041-1042] prepared a graph 
of the dimensionless groups now known as є, R and Nty for 
the 1-2 exchanger. Such a graph had the added versatility 
of simplifying the calculation of performance in a given 
exchanger when operating at conditions different for those 
for which it was designed. Kays and London [Ref. 2: 
pp. 63-74] prepared similar graphs and tables of 
= = f(R, Ntu, flow arrangement) for the 1-2 exchanger and 
for several cases of crossflow and periodic flow. 

The foregoing describes the early history of the search 
for the so-called Logarithmic Mean Temperature Difference 
Correction Factor, F, with regard to heat exchangers having 
an even number of tube passes. It is a fact, however, that 


certain space economies could be realized from exchangers 
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having an odd number of tube passes so that the tube side 
fluid could enter and leave the exchanger at opposite ends 


of the exchanger (see Figure 1.2). 


B. FISCHER'S WORK 

With the foregoing in mind, an extensive search has been 
conducted to obtain e - Nq, data for the so called "1-3" and 
"1-5" exchangers. This search has uncovered a single work, 
that of Fischer [Ref. 11: pp. 377-383], which summarizes 
the historical development covered here and contains only a 
small section on the 1-3 exchanger. This work by Fischer 
develops an equation for true mean temperature difference of 
the 1-3 exchanger and casts the results in terms of F rather 
than с. Moreover, the work treats onlv the case where the 
three tube passes are arranged with two in counterflow and 
one in parallel flow (1-3:2C) making no mention of the one 
counterflow and two parallel flow (1-3:2P) case (see Figures 
2.1 and 2.2). In addition, the equation developed to yield 
F must be solved using a trial and error solution. 

The present work is aimed at continuing the Fischer 
investigation for several reasons: 

1. A solution is needed for effectiveness, e, as a 
function of capacity rate (R) and number of transfer 
units (Nty)- 

2. This solution should be in a closed form if at 
all possible so that it will be computationally 


efficient and useful in both the design and analysis 
frameworks. 
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Figure 2.1 1-3:2C Three Tube Passes - Two in Counterflow 
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3. Because valid design data can evolve from a 
polynomial approximation. The search should not be 
abandoned just because a closed form solution does not 
result from an analytical approach. 


4. Data is needed for the (1-3:2P) two parallel-one 
counterflow configuration. 


5. A 1-3 exchanger in a marine (shipboard) application 
may result in a considerable space saving over its 
l-n counterpart with n even. This would be evident 
on the outside of the exchanger where it would be 
immediately noted that the 1-3 exchanger has tube 
side inlet and outlet at opposite ends of the 
exchanger. 

The next section confirms Fischer's result and shows 
that a closed form solution cannot be obtained for the 
effectiveness of the 1-3 exchanger. Sections IV and V 
demonstate how, through numerical analysis assisted by a 


computer, a polynomial solution can be derived that will 


yield the effectiveness to engineering accuracy. 
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III. AN ATTEMPT AT A CLOSED FORM SOLUTION 

A.  EFFECTIVENESS AS A FUNCTION OF CAPACITY RATES AND 

EXCHANGER SIZE 

This section deals with an investigation into the 
effectiveness, e, of a one shell pass and three tube pass 
heat exchanger, whereby £ compares the actual heat transfer 
rate to the thermodynamically limited, maximum possible heat 
transfer rate as would be realized only in a counter flow 
heat exchanger of infinite transfer area. This exchanger 


heat transfer effectiveness is given бу 





_ а _ Са ров, за“ Thot ,out? a Celteold, out” teold,in? 
P CC OMM cw o IINE CENE ZU - t 
Imax min^ hot,in  “cold,in min' hot,in  '"cold,in) 


where Сујп is the smaller of the Cp and Coe magnitudes. 

Thus, € possesses the significance of effectiveness of the 

heat exchanger from a thermodynamic point of view, with the 

magnitude of the effectiveness completely defining the heat 

transfer performance. In general we express e = f(Ntu, R, 

and flow arrangement) and when the flow arrangement is under- 

stood, it is said that є = f(Ntu, R). [Ref. 2: рр. 14-26]. 
The number of heat transfer units Ntuų is a nondimensional 

expression of the "heat transfer size" of the exchanger. 

When Nc, is small the exchanger effectiveness is low, and 


when Nr, is large, e approaches the limit imposed by the flow 
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arrangement and thermodynamic conditions asymptotically. 


From inspection of the definition of Ng, 








А 
ТЕЛІ С 
Neu ше. Qu | шаг 
min min $ 


it is clear that the overall conductance and transfer 
area affect the costs of attaining a high value for Мец, 
ergo high є. The capacity rate ratio, R, as defined by 
(299 
в = Jin 


= — 


max 


is simply the ratio of mass flow rate times specific heat 
capacity for the two streams. These can be considered as 
flow stream thermal-capacity rates, i.e., energy storage 
rate in the stream per unit of temperature change. (Нег. 2: 
pp. 14-26] 

The attempt taken in this thesis to develop a closed 
form solution has used the basic fundamentals of heat 
transfer as well as those indicated above. A closed form 
solution for e was sought for both 1-3 exchangers with one 
having two out of three tube passes in parallel flow and the 
other having two out of three tube passes in counterflow. 


The analytical approach taken, and demonstrated in this 


section, is for two out of three tube passes in counterflow. 
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B. ANALYTICAL DEVELOPMENT 

The derivation for the effectiveness, є, of the 1-3 
exchanger as a function of the capacity rate ratio, R, and 
number of transfer units, Nty, depends on several 
assumptions. 


(1) The overall coefficient of heat transfer, U, does 
not vary within the exchanger. 


(2) The specific heat of both hot side and cold side 
fluids does not vary. 


(3) Each fluid is thoroughly mixed, that is, the 
temperature of both hot and cold side fluids is 
uniform over any cross section. 


(4) Steady flow conditions are maintained. 


(5) Heat losses to or from the environment are 
negligible. 


(6) No change of phase takes place; all heat 
transferred is sensible heat. 


(7) There is equal heat transfer surface in each pass. 

The configuration is shown in Figure 3.1 where the three 
tube passes are designated with subscripts a, b and c. The 
temperature of the hot (shell side) fluid is indicated by 
upper case letters. For the cold (tube side) fluid, lower 
case letters are used. The subscript 1 always refers to the 
fluid inlet and the subscript 2 always refers to the fluid 
outlet. 

With Мур and Cph designating mass flow (kg/sec) and 
Specific heat (Joules/kg °K) of hot fluid entering at 171 айы 


leaving at To we define a capacity rate for the hot Side 


Ch = VnCpn 
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Figure 3.1 Two Counter One Parallel Configuration for 
Development of Effectiveness Relationship 
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In similar fashion for the cold side (with We and Сре) 


entering at tj and leaving at to, we have 


Се 7 WcCpe 


We then obtain an energy balance for the entire 


exchanger 
сиСТЕ - То) = Celt? TD (m 


Over the right hand side of the exchanger (Figure 3.1) 


ОНО = T) = Colt - te ИИИ (2) 


and a differentiation gives 


Срат = с Саве атр ас) (3) 


Across dx, with a (m2/m), the surface per running meter 


3aL is the total surface in 


of length of pass so that A 
the exchanger, we may write the heat transferred to the 


element dx in each cold pass. 


Cedtp = -Uadx(T - tp) (49) 
(4с) 


Саба жы Шайх с) 


Here it should be observed that due cognizance has been 


taken of the direction of the flow in each cold fluid pass 
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with respect to the positive sense of the length coordinate, 
x, and U is the overall heat transfer coefficient (W/m2-*C). 


With eqs. (4) in eq. (3) 


CQdT = Ua(3T - ба - th - te)dx 


or 
ат 
ТО ЕСЕ) 05) 
where 
Ua 
n = 
ћ Cn 
is a sort of Nc, per unit length for the hot side. 
Now differentiate eq. (5) 
d 2 сх ах ах ах 
Х 
and with eqs. (4) substituted 
а2т d 
ат 3n, соб п.п; (Т = а + ty = Ба (6) 
ах 
where 
ша ша 
шог, 
с 


where again the resemblence of ne to Куц can be noted. 
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From eq. (2) we obtain 


Cn 
we (lope) = 65 - 5 ағам (2) 
C 


and with eq. (7) put into eq. (6) we obtain 


d^T dT _ _ h _ 2 
—5 Зе a Ge СТ ММ 
ах C 
Or 
агт = 50 сЕ -n_n 2 [(В - 1)T 41, - Е ъ. ра 
4 2 и dx сп СУ с 1 с 12 
Х с 
where 
Ке = Се/Сћ 


is the capacity rate ratio for the cold side. 


Notice that 


and 
R = 1 _ Sh 
h Ra mA 
С с 
a capacity rate ratio for the hot side. Then, algebraic 


adjustment provídes 
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ат ат h и 2 _ 
Бе? - 3n, dz +m eas = ng Ti) (9) 


which is a linear, non-homogeneous, second order 
differential equation with constant coefficients having a 


complementary function 


Т = Се + Се (10) 


where Cy and Co are arbitrary constants and where 





Sn 1-8 1/2 
e h 1 2 h 
орнар у tg ба л 
E 3n, Dh T» "m E - 1/2 
2 зид и 2 Е, 
h 
But 
0с с. 2 
m _ (Ua) : h - als = р ^ne eel 
"OP xe denn DS Ge ZES ^ >> 
n, (С) (Па) с На 
со that 
а. „а, = а $ ah [9 - 48 і и 2 B 
li» DO аа h Ry 
or 
Ор 
01599 = з (or o КУ GLT) 


з} 


where 


х = [9 - 4R (1 - R9]? (12) 
Designate the particular integral as Tp, and by the 
method of undetermined coefficients let T5, - P so that in 
eq. (9) 
Do ї 
Р 2 _ 
т ( В ) P = m T4) 
h 
This makes 
Ё В 
Ж и 2211 h 
Tor 2 ай И га 
h h 
so that 
Е ta 
pi - (13) 
са В, 


The general solution to eq. (9) is the sum of eqs. (10) and 


(13) 
Oa. X Oak Хе чот 
T(x) = Cle нг Ce с оо а (14) 
h 


where the arbitary constants, C4 and Co are evaluated from 
conditions at x = 0 and x = L. At x =O, T(x = 0) = To and 
at x = L, T(x = L) = Ту. When these are inserted, in turn, 
into eq. (14), one obtains a pair of linear algebraic 


equations in the unknowns Cj and Co 
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where Tp, is given by eq. (21277 


It is only a matter of algebra to show that 


aS L 
OMEN Sese 
с = 1 рї 2 рї (15a) 
1 a,b ag L 
e ZT eE 
and 
а 
Bd (То - noD - (T4 - E T 
Dm Са a, L (195) 
1 2 
е LO 


It is easy to see from eq. (1) that 


SO that 


to itx c Te 


Use of this in eq. (13) shows that 


ET 


t. + RA CT, - То) - R T4 


Вы. = : 


pi _ 
аг 
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ОГ 


75 


T. = 
рі I-R- (16) 
indicating two alternative forms for Tpi given by eqs. (13) 
апа (16). 


Insertion of eqs. (13) and (16) in eas. (15) for C, and 


Co will yield after some algebra 


11 - to Та t4, «b 
ime c eeu T 
ME ae 1 (17a) 
1 nd aL 
e - е 
and 
т==г о, T. - t 
2 1 1 
белы. на (Tk? 
с. = (175) 
2 a L а 
е - е 


Equation (14) is an expression for the hot side 
temperature at any location in the exchanger in terms of the 
extreme temperatures, t4, to, Tj and То. 

Next take eq. (5) and set it equal to the derivative 
of eq. (14) noting that C1, Co and Tpi are all known 
constants. 


TN е са ЗЭС 
4х = n, (3T = За зли = ) = ПЗ + ag C 4e е9 


At x = 0, where T = To, tag = tj and tp = te = the 
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ат 


and if we subtract eq. (4a) from eq. (4c) we obtain 


пр, aaa 
a © 


ve t 
a с С 


which can be integrated using C3 as the constant of 


integration. 


and at x = O where tag = ty, and te = the 


ti - tbe = СО 


or 

tbe = t1 - ©З 

In addition at x = L, t = tab and са = 12 50 that 
n, L 
E с 

t ig to = Све 
or 

Qu ab 220 Ee 

3 “nb ah 2 
е 


This gives a relationship between tgp and tpe 


N 


= e 2 с 
t = t NS Соје (20) 


9 


where Ne = псЬ сап be considered as the total number of 
transfer units for the cold side. | 

Return now to eq. (18) and look at the conditions at 
x = L where ta = їр = Тар, te = tg and T = T,. These 
conditions in eq. (18) give 


eL On L 
Cie + An Coe 


по (ЗТ, - 21 у - to) = ај 


where again we remember that C4 and Co are known constants. 


Solving for tap 


L L 
_ 1 “1 - 
and with this in eq. (20) 
N a, L L 
с 1 2 Е 


Then with eq. (21) in ea. (19) 


N 


с 
то с ге = t5)] 


0161 202 = Зпу I(T. 


N L І, 
-е ECC * „ее “ ) (228 


Equation (22) confirms Fischer's result [Ref. 12: 
ор. 377-383] and at this point in his development he 


branches off to seek an expression for the Logarithmic Mean 
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Temperature Difference Correction Factor, F. The attention 
here is focused on Е = f(R, Ney) and the balance of this 
section continues in this vein. 

Look at the ajCy + agCo term in eq. (22). Use of eq. 


(11) allows the representation of а1С1 + a9Co 





Пр, Пү 
то (3 t 325, + бе (3 - С. (23a) 
or 
Әп, у 
5— (C, + Co) + -а (С. - С.) (23b) 
Then from eqs. (17) 
a, L L 
(То = tare че е 2 ] 
CES NICA o o 
1 2 а. Г. aL 
(1 - В Ме ЗЕ | 
or 
Ta~ f 
ME хо l (24) 
1 - Ry 
Moreover 
L p 
ГТ. Е.) (Т % ye | - е 2 
Ё 1 2 2 1 
С, ы Со єана а аш ыа зт а Т. Oy L (25) 
(1 - В Ме даде | 
Now let 
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За 


БЭ? 


ш = (26а) 
апа 
2 = нэг (26b) 
2 
SO that 
a, L QoL 
+ 
" [o E Д a ща 2). _ ес 72) 
m „1. 21, E е 721, 
и ЗЕЕ а 2220) 
ОГ 
а: 16 Gol 
c Lc ИЕ (27) 


Moreover, it is easy to see that 


e + е - 2e E cosh zL (28) 


If eqs. (24) through (28) are collected and put into the 


expression of eq. (23b), the result is 


3n, му 
с С. + ч2Со = (C4 + Со) t om (С. - Co) (23b) 
ul, 
To cnp Z(T. —— t.) = Пе esch zL 
2 1 1 2 2 1 
u el H a l Ј 
17 x В, 22е sinh zL 
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ОГ 


пато бо 
ар Са + а Со < 8(-4-т- 50 8 
T4, - t 


2 


ші сер 7 = эг) СООН 21] (29а) 
Ћ 


T. - t 
21 С-г-г-н 067 
and this could also be written as 


Tus rm rb 


DOCERE - 
b Вана уа овен ZL (29b) 
1 В, 


The next step is to reduce the right hand side of eq. 


(22). Use of eq. (11) permits the representation 


2-1: ао Зп An on An L 
il 2 h h a, T. h h “1 
а „Се + «Се = = EE ТЕ (=== - —5-)€5e 
or with eqs. (26) inserted 
а L а L а Т, Г, L L 
1 эс 1 52 pil 2 
а Се t aS Ce - «(Се + Coe E “(Се = Coe ) 
(30) 
But by eqs. (11) and (17) 
со Is (EST (С 9) 15 
а. Г оре Пра а = << о (пе на" е 2 Сг) [е 1772 522 
1 2 12222 21 
с.е + Cae = 
1 2 a, L aL 
(1-R)(e зе") 


ОГ 


Сүе T Coe = ie (327) 
аа QoL 
Сүе - Coe may also be simplified. Again using 


eqs. (11), (17a) and (17b) 
L L ( a t as )L ( азор 18 
Т, кое 1 + е 2 )-(T57t, Te 1 2 + е 1 7 ] 


L Е 
(1 – пре | Ен 


Но 


The exponential term at the far right in the numerator 


is really quite simple. From eq. (11) 


and by eq. (26а) а] + а2 = Зпн = 2%. Thus with the 
combination of exponentials in the numerator and the 


denominator given by eas. (27) and (28) we find that 


20 aL 2(Т, - взе? cosh zL - (Ta - t. )2e 
1 Да 2 1 а) е 2 1 
Ио Е al, 
2(1 – Ве“ sinh zL 
or 
a, L GL T t To Xu 
1 gi Ті 2 o D ie. 8 
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Now with eqs. (31) and (32) put into eq. (30) we 


obtain 
а, ao L Ш ғ 
1 ob 1 2 
а Се T аоСое = шш 
DM = Е 
+ >) coth 21, - ки с. csch zL] 
"n h 
or 
a,C,e + аг Со? 5 к=н) | + Z COCH ZL] 


ore F 


Ян ое )° csch zL (33) 
With eqs. (29b) and (33) inserted into eq. (22) we 
obtain 
Ta = t шэг ae 
P | шо = ео? zL] + S йн csch zL = 
h h 
Но 
За Г СТ, = 51) + е (Т. = Со) = 
=N тит Deut 
e Дас сота 12) = ЕС, zLi 
h h 
(34) 
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We wish to develop an expression for the exchanger 


effectiveness so we designate the hot side effectiveness as 


Tom 
с. ЩЕ 1 2 (35 
h Ti - t4 ) 


and we begin by simplifying eq. (34) by dividing throughout 


by (To - t1)/(1 - Bh) to obtain 


12-21 N Т. e 
-V 
(w = 2- -Coth 4) t Ze ) csch Zit 2 = = КЕ + е орос р 
т. = Е Та то б 
2 ! 2 1 
-Х TEM sb 
= сео сои) маан = ze" csch Z) (22 
2 1 
where 
W = wh (37a) 
Z = 21 (328 
апа 
R = 3np(1 - Rp) (376D 


We can then let 


Ti 7 to 


0 = т 
To 7 ty 


and get eq. (36) to look like 
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(0 -28СО0ТЛ 2 (ге! csch 2) ф5 


ЈА -N. (И-М) 
В + (е “в)ф - Ге ©°( * coth 2)16$ - (ze C esch Z) 


It is now a matter of algebra to solve for ф 


02-22 
_ В + 2е сера (о —- 2 СОН Я) о 
Е и ры ото 5. 
ze" csch Z - Ве + е СО t Z Јуве A) 
The next step is to represent $ as a function of єн. 
This is done with some algebraic gymnastics as follows: 
ака 2. Шле и шинж” 
То - t То -- Т, + Ti - С 2 d - ка 
1 ЧЕ t 
1 
ог 
ти ес 
(Т, - t) - е) 
Могеоуег 
T4 - te Т, ты t 7 to 


UT crema mcm тул т 
(Ty са) (1 3% (71 - "ода е) 
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с Про (72 Дині 
Ty 7.73 H7 T 


1 - е 


1 - 


But Ву = (Солета 00101122)  -о2эПї 


1 = ey Ry, 
"Лтд 
Or 
nep E 
бо SUE E 


where ф 15 given by eq. (38) 

The neatness of the form of eq. (39) is deceptive 
because, unfortunately, it cannot be used to determine a 
unique value for єр. The reason for this can be found in an 
inspection of eq. (38) which provides the value of 9$ which 
is used in eq. (39). 

Notice in eq. (38) that 7, W and Ne are all functions of 
the product aL. On the other band, w, z and R are functions 
of a only. Thus, it is impossible to vary a and L indepen- 
dently and still achieve a unique solution. 

For example, suppose a = 50 m@/m and L = 5 м so that 
aL = 250. A value of ф may be obtained from eq. (38) using 
these values. However, if a = 100 and L = 2.5 so that aL is 
still equal to 250, an entirely different value of ф 15 


obtained because a = 100 rather than 50. 
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One should resist the temptation to multiply numerator 
and denominator of eq. (38) by L thereby creating a 
situation where only Z, W and Ne appear along with a new 
R' = RL. Such a procedure is doomed to failure because in 
dealing with an equation derived from n equations in ntl 
unknowns, one cannot create the п+1 В equation by 
multiplying one of the n equations by a constant. This 
makes the ntit equation so obtained linearly dependent on 
one of the original n equations and the entire set becomes 


linearly dependent. 


This section represents an attempt to obtain є = f(R,Nty) 


and the attempt has not been successful. It is now time to 


turn to the computer and this will be done in Section IV. 


47 


IV. NUMERICAL AND COMPUTER ANALYSIS 


With a closed form solution for e as a function of R 
and N+, not attainable as indicated in Section III, it 
becomes apparent that an alternative method is needed to 
determine the effectiveness for the 1-3:2C and 1-3:2P heat 
exchangers. Kern and Kraus [Ref. 12: pp. 306-360] describe 
a computer code for a thermal analyzer. This code (program) 
makes use of node equations generated by finite differences 
and it employs a Cholesky LU decomposition scheme. 

The Cholesky decomposition, as explained by Stewart 
[Ref. 13: pp. 134-144], is best used when decomposition in 
the presence of positive definite matrices is requested. 
This is the case at hand where one tries to solve 
numerically for the temperatures that lead to the 


effectiveness of the 1-3:2C and 1-3:2P heat exchangers. 


А. THERMAL ANALYZER TVSSI 

The computer program employed is Program TVSSI (Appendix 
A) which is an adaptation of the thermal analyzer program 
called TVSS2 and listed by Kern and Kraus [Ref. 12: pp. 
306-360]. The adaptation consisted of changing the program 
so that it could perform the computations in the SI system 
of units and be receptive to the use of a Specially created 


input file. Also it should be noted that TVSS2 was written 
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to be used in conjunction with the Honeywell H-1800 computer 
system. Therefore, it had to be modified to run on the 
Naval Postgraduate School's IBM 3033 AP system. 

The program itself is a non-linear equation solver that 
determines the temperatures at a prescribed number of node- 
points or nodes from a set of node equations in almost any 
framework (i.e., network analysis, field plotting, or fluid 
flow distribution). It has certain features that make it 
primarily an equation solver for thermal analysis. These 
features include: 

1. an ability to linearize radiation terms. 


2. an ability to allow any of the coefficients in the 
node equations to vary with temperature. 


3. an ability to provide constant heat input and heat 
input as a function of temperature at any node. 


4. an ability to consider other modes of heat transfer 
that are non-linear such as boiling and natural 
convection. 

As stated earlier, the program utilizes the Cholesky 
decomposition scheme and, because of the linearization of 
the radiation terms (a feature of the program that is used 
even though radiation does not appear in this е - Ntu 
Study), the program is iterative. 

Cholesky's decomposition consists of finding a lower 


triangular matrix [L] which is capable of reducing the 


original system of equations. 


2 | 7-0) (40) 


49 


ОГ 


ГЕТІДІ Со OI (41) 


to the unit triangular form 


АТА — ПВО (42) 


SO that the sought after elements of the column vector TT] 
can be obtained by backward substitution. 

Suppose, for example, that [K] is 3x3 and assume that 
the system [K][T] = [Q] has heen reduced to the form 
ГАЈ[Т] - [В] = [0]. In this event a premultiplication by 


[L] will return the system to its original form, that is 


[L]CTATTT ] PB ГУИ 1 


This implies that 


[L][A] [K] (459 
and 


[LJTR] 


[0] (44) 


These equations allow the determination of [I], [A], and 
[B] in a very simple manner and the matrices are uniquely 
determined because [K] and fO] are known, or, at least are 
known after each iteration because the elements of [K] are 


linearized. For a 3x3 system 
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[K,Q] [L] [А,В] 


К11 К12 Кіз ај 111 0 0 В! 
Кој Koo Коз 92| =| 127 122 0 ин! a23 bo 
k31 Кзо k33 аз 2 о 1541 10 0 1 b3 
one may obtain the following for the elements of [L], [A], 
and [B]. 
k11 ? (1)141 + (0)142 + (09143 = lil (45) 


which shows that the first column of [L] is identical to the 


first column of [Kl]. 
шэн. їг 200022-0002-34/27-11181|) = Кітаз і 6) 


which shows that the first row of [А] is equal to the first 


row of [К] divided by [k11] and then 


КР oma el) “7ЄЯВ Бо р2529Кво. -7То18:19 

Ын 2 121:1:21!22222 222: = (КЗ 1о1а183/125 
ЕИО ЕСЕ ЮКЕ 13 еве” Палацо 

Ср о Шо о 2022-2002 -41210121192 


In the foregoing manner the elements of ||, [A] and [В] 
are obtained successively in terms of previously determined 
elements in a progression that goes horizontally from 122 


on. Thus the general relationship are seen to be 


S 


Ли! 


14-37 Ж 15) (47) 
г=] 
with 
MSS (48) 
and 
i-1 
т = 
aig on, Ses 11565! СО) 
qr r=1 
with 
к.. 
DES (50) 
ij Kaa 


Moreover, it is observed that if [K] is symmetrical which it 


must be in our coupled set of equations (kij - Кјј) then 
1114 з ajilii (1; j = 1,2,3 ma mn (51) 


The modification for computation in the SI system was 
quick and simple. It involved changing a numeral in two 
places (460 to 273) and some format statements (°F to °C and 
Btu/hr to Watts). 

The conversion of TVSS2 to TVSSI for running on the 


IBM 3033 AP system (which uses the FORTVS compiler) required 
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a modification to the Fortran program language used in TVSS2 
in order to compile the FORTVS (basically international 


Fortran 77) system used on the IBM 3033. 


В. INITIAL MODELING 

Initially, the model was designed to find the tempera- 
tures that would allow computation of a single effectiveness 
value after a detailed set of capacity rate, coefficients 
and surface data were entered into the program. From the 
scope of the problem it was realized that multiple runs of 
the thermal analyzer program (TVSSI) would be needed. It 
therefore became necessary to develop a program that given 
Ch, Ce, U, A, T1, and t1, an input file would be created for 
use by the modified version of the thermal analyzer 
(27551). 

The first step taken was to develop a program to create 
an input file for TVSSI that would yield the effectiveness 
for a 1-4 heat exchanger which could be compared to the 
existing analytical solution for the effectiveness of a 1-4 
exchanger. With this accomplished and confidence 
established, similar programs for the 1-3:2C and 1-3:2P 
exchangers could be developed. This program was called 
NTU14 (See Appendix B) and the following parameters were 
used for all runs. 

1, 250 nodes were used. 


2. The initial temperature for the computer to begin the 
iterative process was set at 200°C. 


ES 


3. An eventual accuracy of .05 between the final and next 
to last iterations was used. 


abs A radiation coefficient convergence factor of 0.66667 
between iterations was used. 


5. The maximum number of iterations that the computer 
was allowed to perform was set at 12. 


6. A damping factor of .8 was set as an initial damping 
based on the number of non-linear terms in all of the 
node equations. 

When the values of Cy, Со, О, А, Ту and tj are set, Ми 
and R are then compiled and an input file for TVSSI was 
generated. In this file all node equations and internode 
conductance values were determined. This program determined 
and specified the nodes that interact with each other and 
the methods by which the interaction takes place such as 
conduction, forced convection, and fluid flow. 

The program, NTU14, makes use of the fact that each term 


in a node equation shows three things. The first is the 


node that is coupled for heat flow with the node in 


question. The second is the method of heat flow between the 
nodes. In this case forced convection and fluid flow are 
used. Finally, the node equation shows the magnitude of the 


internode heat flow. Here all the pieces of information are 
collected and presented for use by TVSSI as an input file 
with all items in the proper format. 

A comparison of the effectiveness for the 1-4 exchanger 
developed by the computer to that of using the closed form 


analytical solution for effectiveness developed by Kraus 
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and Kern [Ref. 8] as shown by equation (52) 


o M o M a (<) 


щори 10 й Ан е Mate E аталы | 


4 


was then undertaken. The results of this comparison showed 
that over the entire range of R from .01 to 1.0 for varying 
values of Ney from O to 3.25 less than a 0.5% difference was 
ever realized. A small sample of these results are provided 
in Table 1. The conclusion to be drawn here, is that the 
methodology used to develop the computer program NTU14 for 
input to TVSSI for finding effectiveness was sound and could 
then be used in the development of the 1-3:2C and 1-3:2P 


exchanger methodology. 


C. DEVELOPED MODELS FOR 1-3:2C AND 1-3:2P HEAT EXCHANGERS 
The same technique used in developing the program NTU14 
was used to generate computer programs NTU32C and NTU32P. 
These are listed in Appendices C and D. The departure for 
each of these programs from the NTU14 program is in the 
number of nodes; they are based on 200 node models as shown 
in Figures 4.1 and 4.2. An example of the output file 
generated from one of these programs is found at Appendix F. 
іс these values shown in Appendix F that are used hy the 
thermal analyzer to determine the temperatures To and to 
for the specific set of given initial parameters Су, Сб, К, 


ST, and ti. 
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ТАВЬЕ 1 
COMPUTER TO ANALYTICAL COMPARISON 


РОК В = .5 

Е нэ ANALYTICAL | |. COMPUTER | 006 27-74 

cu RESULTS RESULTS DIFFERENCE 
0.05 ‚0482 0482 0.00 
0.25 .2094 .2090 0.20 
0.50 .3569 .3559 0.28 
0.75 .4628 .4612 0.35 
1.00 .5398 .5377 0.39 
1.25 .5963 .5940 0.39 
1.50 .6379 .6354 0.40 
2.00 ‚6915 . 6886 0.42 
2.50 .7206 Ул 0.40 
3.00 .7360 ‚7333 0.37 
3.25 .7406 ‚7381 0.34 





D. SCOPE OF COMPUTER ANALYSIS 

At this point it is possible to let the computer solve 
for temperature values that yield a value for effectiveness 
based on a specific set of initial parameters. However, 
it must be realized that many computer runs are required to 
generate enough data to ensure confidence in the results 
which cover a wide range of capacity rate ratios and Nry 
values. 

To efficiently expedite the computer task, the Multiple 
Virtual System (MVS) with Job Entry Subsystem and Networking 
(JES3) was utilized. The MVS coupled with JES3 is more 
commonly referred to as batch processing. Based on trial 
and error, it was determined that in order to build a solid 
data base, eleven different values for effectiveness were 
needed to best represent a particular value of R. This 
was required over a range of R from R = 0.1 to 1.0 in 
вповетепсє ої .01. In all, 200 curves for both the 1-3:2C 
and 1-3:2P exchangers were needed. This means that 2,200 
unique values of effectiveness needed to be found, plus 
100 values of effectiveness to be used for comparison with 
the 1-4 exchangers. 

To complete this task, TVSSI was slightly modified in 
accordance with the appropriate guidelines of the job 
control language (JCL) needed to run on the batch processing 
system. These modifications are few and were needed only 


at the beginning and end of TVSSI. The modified version 


БИ 


of TVSSI has been called TVCOUNT with changes shown in 
Appendix F. It was TVCOUNT that was then used to activate 
ТУ. 

It also became necessary to modify the three input file 
programs NTU14, NTU32P and NTU32C sueh that they needed to 
be compiled only once. They were then loaded in a library 
file to be used when called by another program. New 
programs utilizing the batch system were written that could 
easily be loaded with the appropriate input data for a 
Specific R value. These, which are referred to as "sister 
programs," were used to go from the library file to TVSSI 
and cause TVSSI to be executed eleven times under TVCOUNT 
covering the desired range of Nty for a specific R value. 
The revised input files called NTU14C, NTUS32CC, NTU32PC and 
their associated "sister execution programs," NTU14L, 
NTU32CL, NTU32PL, are found in Appendices G through L. 

The overall system flow chart of how all of the fore- 
going is accomplished is found in Figures 4.3, 4.4 and 4.5. 
It is noted from these figures that TVSSI is referred to as 
TVSSIA through TVSSIV. These are the same programs as TVSSI 
but for bookkeeping purposes by the computer they are 


labeled A through V. 


Е. COMPUTER RESULTS 
Upon completion of all data collection from the computer 


output, plots of effectiveness vs. Nt, for the whole range 
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of R were plotted and are shown in Appendices M and N. 
Thirty-three different plotting programs utilizing the 
"Display Integrated Software System and Plotting Language 
(DISSPLA)," were written to graph the data obtained. 

An example of one of these programs is provided in 
Appendix O. 

Examination of the graphs in Figures 4.6, 4.7, and 4.8 
shows that the 1-3:2C exchangers outperform the 1-3:2P, 1-2 
and 1-4 exchangers. This is true in all cases, and, because 
of this, only a sample of the data was chosen to be shown in 
these figures. Furthermore, it is noted that at higher Nr, 
values, the effectiveness of the 1-3:2C exchanger is better 
than all of the others considered, while at higher capacity 
rate ratios, the effectiveness of 1-3:2C exchanger even 
begins to outperform the others at lower Nt, values. This 
increase in performance is easily understood because it has 
been proven by Kern [Ref. 1: pp. 139-137] and others that 
greater temperature differences result when process streams 
are in counterflow than parallel flow. Thus, when there is 
a combination of the two phenomena (counterflow and parallel 
flow) occurring, then it becomes apparent that the extra 
counterflow pass must increase the exchanger's overall 
performance. 

As shown in the graphs in Appendices M and N, effective- 
ness increases as both R and Ng, increase. These curves can 


be used to give a graphical approximation of effectiveness 
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if that is all that is required. From the empirical data 
used to develop these curves, further investigation into the 
development of equations for these curves which may be used 
to determine an exact value of effectiveness when R and Мет 


are known, is undertaken in Section V. 


60 


EIE (515141 ` 
(1 55у4) а14 алоо 4 
«о 
Bae вы 


$ (г $5\а) апла ачоов 


«(ttt 


«(€ Ssva) aima a1024 


$ aina LOH È 





13G3OWN 3GON Э26-1 


Ен п 


Nodal Arrangement for the 1-3 


Figure 4.1 


ow 





В (1 55у4) апза 0709 $È 





4 (е взуа) аплч 9109 4 


b (€ 65у4) апла ачооф 


«Ста CLE: 


È ama Lon $ 


забори заом ас'є-і 


2P Exchanger 


Nodal Arrangement for the 1-3 


Figure 4.2 


62 


II юэ of TR parameters 

given Си; Се, О, А, Ту апа tj) 
make ир t input” data to be 

solved for 


"ed 
-Generates input file 
NTU14 [|n TU14BC under Batch Processing 
-Stores file in a Library 


Processing 


-Uses Batch Processing 


-Calls input data 
-Calls Library input file 
NTU14BL -Calls TVSSI 11 different 
times and redesignates 
TVSSI as TVSSIA thru TVSSIK 


-Modified thermal analyzer 
VCOUNT -Uses Batch Processing 
-Executes TVSSIA thru TVSSIK 


-11 different corresponding 
values for effectiveness 
based on initial input data 
is listed here 


Figure 4.3 Computer Systems Flow Chart for 1-4 Exchanger 
Model 
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ll sets of initial parameters 

given (Мец Ch» Cc, U, А, Ті апа tip 
make up the input data to be 

solved for 


Modified 


for -Generates input file 
МТОЗ2С а NTU32CC under Batch Processing 
Processing -Stores file in a Library 


-Uses Batch Processing 


-Calls input data 
-Calls Library input file 
NTUS2CL -Calls TVSSI 11 different 
times and redesignates 
TVSSI as TVSSIA thru TVSSIK 


-Modified thermal analyzer 
T -Uses Batch Processing 
VCOUN -Executes TVSSIA thru TVSSIK 
-11 different corresponding 
values for effectiveness 
based on initial input data 
is listed here 


Figure 4.4 Computer systems Flow Chart for 1-3:2C Exchanger 
Model 
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ll sets of initial parameters 

given (Nt, Ср» Сс, О, А, Ту апа t1) 
make up the input data to be 

solved for 


-Generates input file 
NTU32PC under Batch Processing 
-Stores file in a Library 


-Uses Batch Processing 
-Calls input data 


-Calls Library input file 
TUS32PL -Са115 TVSSI 11 different 
times and designates 
TVSSI as 17550. thru TVSSTV 


-Modified thermal analyzer 
VCOUNT -Uses Batch Processing 
С -Executes TVSSIA thru TVSSIK 


-11 different corresponding 
values for effectiveness 
based on initial input data 
is listed here 





Figure 4.5 Computer Systems Flow Chart for 1-3:2P Exchanger 
Model 
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Figure 4.6 Comparison of Analytical (1-2 and 1-4) to Computer 
Results (1-3:2С апа 1-3:2Р) а к=к? 
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Figure 4.7 Comparison of Analytical (1-2 and 1-4) to Computer 
Results (1-3:2C and 1-3:2P) at R - .5 
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Figure 4.8 Comparison of Analytical (1-2 and 1-4) to Computer 
Results (l-3:2C and 1l-3:2P) at R = 1.0 
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V.  POLYNOMIAL REGRESSION 


A. DEVELOPMENT OF POLYNOMIAL EQUATIONS 

The empirical data obtained for the two 1-3 heat 
exchangers was designed to cover an extensive range of R 
pes varying from 0.01 to 1.90 in increments of 0.01. As 
discussed earlier in Section IV, the data obtained at a 
specific value of R is the computer evaluated result of the 
ENUEetiveness, for an associated Nt, value. With this 
accomplished, it then hecomes possible to graph separate 
curves for each of the different R values as shown in 
Appendices M and N. Through a polynomial regression 
technique, as discussed in this section, it is also possible 
to develop implicit equations for the curves with е Ї(Мь, 
R). It is also apparent from an inspection of the graphical 
representation of the empirical data in Anpendices M and N, 
that the curves conform to a high degree polynomial. 
However, further analytical investigation is needed to 
ascertain the exact order of the polynomial terms. This 
investigation will not only lead to the order of the 
nolynomial, but to the specific equation for each curve. 

By use of polynomial regression, the least-squares 
method can be readily extended to best fit the data to the 


mth-degree for the polynomial 


Na 


у = Ас + Ајх + Аох2 г... А (53 
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with the error defined by 
ei = 05 – уз = 0, – Абд - Арх з ДАохі? - ... = АВ 


where Dj represents the empirical data value corresponding 
to xj, xj being free of error. 
The objective is to minimize the sum of the squares of 


the residuals, Sr, 


a 2 n 2 m 2 
524 - E є = 5 (Di = Ay = АХ. = А.х, ee А, ) (54) 
Because at a minimum, the partial derivatives 98г/ ЗА0, 
98p/93A1 ... 95у/9Ар vanish, after taking the derivative of 


Sp with respect to each of the coefficients of the роју= 


nomial, it can be seen that 


Then by dividing by -2 and rearranging we obtain 
АОМ + Ау ХУ xq + A1 x42 ишийг тш xao v» Da 
Ag E Xj + Ay E xi^ + До Үү xi? + ss. А 2 x,m*l zo NAT 


Ag I xj? * Aj ZI xj? * Ao Еэ хү3 * ... * Ag I x,"*2 - у; х,2р, 


Ag E xj + А] хх: МТ + Ао г х #2 +... + Ад Ех 20 = у хур, 


where all summations are from i=1 through n. All of the 
foregoing т+1 equations are linear and have m+l1 unknowns: 
Ао, Ај, Ао, ... Аһ. The coefficients of the unknowns can be 
calculated directly from the observed data. Thus, the pro- 
blem of determining a least-squares polynomial of degree m 


is equivalent to solving a system of m+l simultaneous linear 


equations. Putting the equations in matrix form yields 

N L Xi L хі? L xi? и Du 

У хі L x42 2 xi? L Хас 22” ха 1+1 У ху) 
г хі? TE ;х.4 го хі . . . є хі ПФ РА) -|у хі2Пі 
КІШ o xnl у; х,т"2 5 ху 093 , . , g x,2m ы халі 





[Ref. 15: pp. 302-309 and Ref. 16: 468-474]. 
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From this point on, one finds that it is best to use a 
computer to assist in solving the simultaneous equations and 
this will also help alleviate any ill-conditioning that may 
otherwise occur. An existing "curvefit" program available 
through NON-IMSL [Ref. 16] and found in Appendix P was 
used although some modifications were made to the original 
program to best accommodate the goals of this work. 

To determine the order of polynomial that should 
eventually be used, one increases the degree of the 
approximating polynomial as long as there is a statistically 
significant decrease in the variance c? , which is computed 


by 


ал уш эы ш шы (55) 


In otherwords, the selection of the optimum degree poly- 
nomial is contingent upon a decreasing variance and once the 
variance begins to increase, the degree of the polynomial 
becomes too high. For all cases, it was found that the 
€ - Ney developed curves are of the 5th order. 

As shown in Figures 5.1 and 5.2 the computed values of 
effectiveness vs. Ny for R = 0.1, 0.5 апа 1.0 for both 
flow arrangements, (1-3:2P) and (1-3:2C), have been graphed 
and fitted by a 5th degree polynomial. Because all computed 
values for effectiveness follow a predictable trend, 


only a sample of the data covering the whole range of 


122 


values of R have been shown. It is clear that the graphic 
interpretation strongly backs what is known analytically 
from the polynomial regression technique. Where the 


меша слопоћјр Тог Е = (Мец, ER) is found explicitly from 
є = As5Nty? + AgNtu* + A3ZNtu? + AgNtu2 + A1Ntu + Ao (56) 


while the corresponding coefficients A5, Ад, Аҙ, Ао, А), апа 
Ag relating to a specific value of R are found in Tables 2 
Па тог the (1-3:2P) and (1-3:2C) configurations. Ап 
example of how to use this equation in a heat exchanger 


problem now follows. 


B. NUMERICAL EXAMPLE 

Consider a heat exchanger containing 400 m2 of (A = 
400 m= ) of surface and operating with an overall heat 
transfer coefficient of 80 W/m2°C. Cold fluid at a capacity 
rate of 10,000 W/°C enters the exchanger at 60°C. Hot fluid 
at a capacity rate of 20,000 W/°C enters the exchanger at 
200°C. 

1. Find 

The effectiveness (£) and compute the hot and cold 

fluid outlet temperatures for the (1-3:2C) shell tube pass 
arrangement. 

2. Assumptions 


1) Negligible heat loss to surroundings and kinetic and 
potential energy changes. 


2 


2) Constant thermal and fluid properties for both 
fluids. 


3. Analysis 
Here Ce = 10,000 W/°C and Cy = 20,000 W/°C this 


makes R = Cmin/Cmax = Ce/Ch 7» (T1 - T2)/(to - t1) 


10,000/20,000 = 0.5. 


and Ney ЏА /Се 


80 (400)/(10,000) 


! 
Од 
N 


First, go to Table 2 (page 81) for the (1-3:2C) arrangement 


with R = 0.5 and find the coefficients 
Ag = 0.1294 x 1072 
41 = 0.98120 
Ао = -0.66161 
Аз = 0.27938 
Ад = -0.66456 х 10-1 


As = 0.66069 x 1072 


Then apply equation (56) for Nyy = 3.2 
c 7 AsNty? + AgNtu4 + AgNtut + AgNtu? + AiNtu + Ag (56) 


с + 0.769 
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Because Ce < Ch 


and with T4 - t4 - 200 - 60 = 140°C 


! 


TOSE ti БОСОВ Ug 


0.769 (140) 


! 


ОМАНА 


Finally, the outlet cold fluid temperature to is 


to = 107.7 + +1 
= 107.7 + 60 


= 167.7°C 


and the fluid temperature, То is easily found 


Ат 
2 
Е = = 0.5 
EE 
п (Е) 
ЫГ 
= 145.1°С 
4. Observations 


The primary observation made here is that by using 
the 5th order polynomial equation (56) with the appropriate 


coefficients found in Table 2 or 3, an accurate value for 


T5 


effectiveness can be found thus allowing one to solve for 
many more unknown values or prarameters of the heat 
exchanger (i.e., hot and cold outlet temperatures). 

The other observation that is to be made is that 
when comparing the value for effectiveness computed here 
against the value for a 1-2 or 1-4 exchanger (0.745 and 
0.740 respectively) under the same conditions, one finds 
that there is a significant difference in exchanger 
performance as a function of odd or even tube passes and 
that the 1-3:2C arrangement has a higher effectiveness than 
either the 1-2 or 1-4 arrangement. From inspection of the 
curves for the 1-3:2P exchanger at Figure N-1 or N-6 witn 
R = 0.0 and Ny = 3.2 an approximate value of e = ./15 TS 
Obtained. It is clear that this is also less than that of 
1-3:2C arrangement. Therefore, it is evident that the 
1-3:2C exchanger out-performs not only tne 1-2 and 1-4 
arrangement but also its counterpart the 1-3:2P exchanger by 


3.1%, 3.8% and 7.0% respectfully. 
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2 OUT OF 3 PASSES IN COUNTER FLOW 
DATA POINTS FIT BY 5TH ORDER POLYNOMIAL 


EFFECTIVENESS VS. NTU 


Figure 5.1 
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1-3:2C Data Fit by a 5th Order Polynomial 


22 
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2 OUT OF 3 PASSES IN PARALLEL FLOW 
DATA POINTS FIT BY 6TH ORDER POLYNOMIAL 


EFFECTIVENESS VS. NTU 


Figure 5.2 
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1-3:2P Data Fit by a 5th Order Polynomial 
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VI. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 
Effectiveness values, though not analytically derived, 
can be determined for the 1-3:2C and 1-3:2P heat exchangers 
by utilizing a 5th order polynomial approximation. 
sufficient data now exists to cover a complete range of 
capacity rate ratios for values of Ny, from 0.0 to 3.295. 
From the knowledge of a particular three tube pass heat 
exchanger arrangement and dimensions, fluid flow rates and 
temperature extremes, Ny, and R values may be computed. 
Then the effectiveness may be determined using the 
appropriate coefficients from the tables provided herein. 
After the effectiveness is obtained, there are two choices. 
1) From a knowledge of qmax (see Section III), it is a 
Simple matter to determine the actual heat transfer 
rate from the postulation that 4 є єдпах: 


2) Both fluid outlet temperatures may be computed since 
а = weCpn(t2 - (1). 


It is also apparent from the work done here that the 
1-3:2C exchanger effectiveness outperforms that of the 1-2, 
1-4 and its counter part, the 1-3:2P exchanger, as Nty 
increases. This is true for any and all values of R. 
Therefore, because it is possible to determine the effec- 
tiveness of a 1-3 exchanger which has a higher effectiveness 


(the 1-3:2C arrangement) than that of the 1-2, 1-4 and 


ОЁ 


1-3:2P exchanger particularly at high Nr, levels, the 1-3 
exchanger can now be given full consideration in heat 
exchanger design. This would be especially helpful where 


three tube passes could alleviate a configuration problem. 


B. RECOMMENDATIONS 
The following recommendations are provided for possible 
follow-on projects of a similar nature. 


e Continue study of the є - Куц method for the 1-5 heat 
exchangers. 


e Investigate the possibility of using a linear 
approximation from the 1-3 data to find a periodicity 
of R values at which to develop a 1-5 data base. 


e Develop interactive software to be used on a 
microcomputer that contains both the 1-3 and 1-5 
polynomial approximation coefficients. The procedure of 
entering the values for Му, R, and the type of heat 
exchanger when asked to do so in a menu-driven fashion 
so that effectiveness values can be readily obtained 
is self-evident. 


e After results of 1-5 exchanger analyses are complete, 
investigate the effectiveness for both the 1-3 and 1-5 
heat exchangers experimentally to confirm or refute 
the theoretical work done here. 
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SAMPLE OUTPUT FROM NTU14 COMPUTER INPUT ANALYZER PROGRAM 
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APPENDIX F 
MODIFIED SECTIONS OF THERMAL ANALYZER TO RUN ON BATCH SYSTEM 
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